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Abstract: 

By continuation from the hyperbolic limit of the cardioid billiard we show that there 
is an abundance of bifurcations in the family of limagon billiards. The statistics of these 
bifurcation shows that the size of the stable intervals decreases with approximately the 
same rate as their number increases with the period. In particular, we give numerical 
evidence that arbitrarily close to the cardioid there are elliptic islands due to orbits created 
in saddle node bifurcations. This shows explicitly that if in this one parameter family of 
maps ergodicity occurs for more than one parameter the set of these parameter values 
has a complicated structure. 
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1 Introduction 



In Hamiltonian systems ergodicity is a rather rare property ||T| . Most systems are "mixed" , 
in the sense that there are regions of regular motion of positive measure and regions with 
stochastic motion, which are expected to be of positive measure, too; see for a review 
on the coexistence problem. 

In this paper we investigate the possibility of ergodicity in the family of billiards 
introduced by Robnik [Q, which is given by the simplest non-trivial conformal mapping 
of the unit circle; in polar coordinates the boundary of the family of limagon billiards is 
parametrised by 

p(0) = 1 +ecos(0) , 0e[-7r,7r] , (1) 

where e G [0, 1] is the family parameter. This family of billiards has received a lot of 
attention, especially in the context of quantum chaos, see e.g. |§, |, |^, and the references 
therein. Corresponding to £ = is the integrable billiard in the circle, while e = 1 
gives the ergodic billiard inside the cardioid |0, For a small perturbation of the 

integrable case the KAM theorem implies that most invariant tori persist, whereas for 
the neighborhood of the hyperbolic case there is no general result about persistence of 
non-uniformly hyperbolicity. From this point of view one of the most interesting questions 
for this billiard family is whether it is hyperbolic or even ergodic for some e <1. 

For £ G [0, |[ the boundary of the billiards is convex, C°°-smooth and the curvature is 



strictly positive. As a consequence of a general theorem by Lazutkin |rU| these billiards 
have a collection of caustics near the boundary, which correspond to rotational invariant 
circles in the Poincare map from bounce to bounce. Since rotational invariant circles divide 
the Poincare section into invariant regions these billiards are not ergodic. At e = | the 
boundary has zero curvature at = ±7r and therefore the invariant circles corresponding 
to the caustics near the boundary no longer persist [11]. For e > the billiards are 
no longer convex, and the billiard map becomes discontinuous. From these observations 
and from numerical experiments Robnik suggested that for s > \ the billiards might be 
ergodic 0, see also [0, |13|. It turned out, however, that at least close to e > | this is 



not true because still a family of Liapunov stable periodic orbits exists |T^. This orbit 
is involved in a cascade of period doubling bifurcations, where the extrapolated limit is 
e = 0.5582 . . . WM. A first example of a stable periodic orbit above this value was given 



m 



151 where a bifurcation is discussed for e ~ 0.66948. 
One way to prove that the system is not ergodic is to establish the existence of an 
invariant curve. This is possible by using KAM theory in the neighborhood of an elliptic 
fixed point of the map. A sufficient condition for the existence of an invariant curve in 
the neighborhood of an elliptic fixed point is that the eigenvalue of the linearised map at 
the fixed point is not a third or forth root of unity and that the twist does not vanish, see 



T6| . In this was rigorously established for a certain orbit of the limagon family. In 
the following we will rely on numerical results in which we cannot easily check the twist 
condition. However, as in one parameter families the twist generically vanishes at points 



T7| , the existence of an elliptic orbit for one given parameter value gives strong numerical 
evidence that the map is not ergodic at this parameter. 

As for e larger than the extrapolated limit of the period-doubling cascade there are 
no obvious stable islands in the Poincare section, it might be conjectured that at this 
point all the orbits that are part of the bifurcation tree that is connected to the integrable 
circle billiard have turned unstable. But even if this was true (for which we have not 
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found any counterexamples), it does not imply ergodicity of the system. The reason 
is that there exist stable orbits which are not connected to the integrable limit. They 
appear in saddle-node bifurcations, in which a stable and unstable orbit are created "out 
of nowhere", in principle at any parameter value. At each saddle-node bifurcation the 
newly created elliptic orbit exists in a certain parameter interval until it turns unstable 
in a period doubling bifurcation. We give numerical evidence that there are saddle-node 
bifurcations arbitrarily close to the ergodic limit at e = 1. Moreover, our numerical 
computations show that the sum of the widths of the stable intervals of elliptic orbits 
born in saddle-node bifurcations has considerable measure in parameter space. For a 
related study of the phenomenon of the appearance of elliptic islands in systems close to 
ergodic billiards see [0 or in the standard map with large parameter see . 

However, we cannot rule out that ergodicity holds for some values e < 1, because the 
stable intervals of these orbits are so small that they do not seem to cover the whole interval 
The method we use is continuation of periodic orbits away from the hyperbolic 
limit. At this limit we know all periodic orbits. For the regular orbits, i.e. those orbits 
which never hit the cusp of the cardioid, a binary symbolic dynamics was proven 

in pT| , p2| , and for the cusp orbits a similar coding with an additional letter can be 
used. The idea to continue periodic orbits from hyperbolic a hyperbolic limit is similar 
to continuation from an antiintegrable limit |2^, which has recently be applied to the 



Henon map ||2^. However for an unbounded system like the Henon map, the map can 



actually be proven to be hyperbolic near its limit [B^, PB[. For a system with bounded 



phase space this is rather unlikely due to the conservative Newhouse phenomenon |27 



Moreover, the cardioid billiard is not uniformly hj^erbolic. There is no positive lower 
bound for the Liapunov exponent of certain orbits. Nevertheless, for the cardioid every 
finite periodic orbit is hyperbolic and therefore can be continued. But due to the non- 
uniform hyperbolicity, the system can become non-hyperbolic when perturbed. 



2 Finding and continuing orbits 

Our approach is to continue periodic orbits of the cardioid billiard (e = 1) to smaller 
values of the parameter e. For the continuation we use the discrete variational principle, 
which is numerically much more stable than iterating the map, and also allows for the 
continuation of singular orbits. The discrete analogue of Hamilton's principle states that 
the action 

= ^L(x„x,+i) , (2) 

i=ti 

is stationary on orbits, where L is the discrete Lagrangian of the system. Performing the 
variation 6W = gives a second order difference equations. Assuming the twist condition 
holds the analogue of the Legendre transformation can be used to rewrite it as an area 
preserving map, see e.g. PB[ for a review. 

For billiards the discrete Lagrangian L{x, x') is just the geometric length between the 
two points on the boundary parameterised by x and x', so that W is the total length of 
the orbit. The Poincare map from bounce to bounce of a smooth billiard is a map of 
a finite cylinder to itself. The global winding number of a periodic orbit is given by ^, 
where n is the number of reflections and m counts how often the orbit goes around the 
cylinder. 
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We use polar coordinates (p, 0) to describe the billiard boundary, so that the polar 
angle replaces x. With ([l|) the discrete Lagrangian becomes 



(3) 



For a doubly periodic Lagrangian periodic orbits of period n are critical points of the 
periodic action 



W^{Xn) = L{xo, Xi) + L{xi,X2) H h L{Xn-l,Xo) 



(4) 



where X„ 



[Xo, Xi, 



n-l 



In general the critical points of are found using 



Newton's method applied to the system of equations VW^ = 0. In ^ it was shown 
that all regular orbits, i.e. those which never hit the cusp, of the cardioid are maxima of 
V^P, such that much more stable numerical methods can be used to find them. 

Let us now briefly recall the definition of the symbolic dynamics for the cardioid 
billiard |20[ which gives a complete description of all possible orbits in the cardioid, 
see for proofs. Consider two consecutive points on the boundary with polar 

angles 0,0' G (— 7r,7r]. If = ivr or 0' = ±7r the symbol C is associated. Otherwise 
the corresponding symbol is A if 0' < and B if 0' > 0. The method described in 



22| requires to find initial conditions consistent with a given code word, for which we 
use the following procedure. Rotate the string until it starts with BA. In step 1 choose 
(00, 0i) = {z, S) with z > S > as the initial pair corresponding to the letter A. Here S is 
a small positive number and z will be determined at the end. In step i take 0j_i from the 
previous step, and construct the pair (0j_i, (0i-i — 7r)/2) for A and (0i_i, (0i-i + 7r)/2) 
for B. In the last step n we obtain by construction the letter B, which is the first letter 
of the word. We have to fulfill the conditions 0„_i < z and z > 6, e.g. by the choice 
z = (max(0„_i, 6) + 7r)/2. Note that all 0j are in [— vr, vr] for arbitrary sequences of A and 
B. This shows that we can choose an admissible initial condition for each code word, i.e. 
there is no intrinsic pruning [p2| . 

If the maximization of W for a given symbol sequence leads to a solution for which 
one or more segments run outside of the billiard, this periodic orbit is not realised by the 
system; the orbit is called "f-pruned" It may also occur that the solution possesses 
a point in the singularity, a situation in which we say the orbit is "s-pruned", or in the 



more general context of ||2^ extrinsically pruned. Also combined types of pruning occur, 
see pO[ for examples of the different types of pruning. 

In order to find all orbits for e < 1 we also have to find finite orbits in the cardioid 
billiard, i.e. cusp orbits, which begin and end with an orbit segment in the cusp, since 
most of them turn into regular orbits for e < 1. This might seem surprising at first; see 
below for a proof. The main difference between these singular orbits and the regular ones 
it that there is no reflection condition to be fulfilled for orbit segments in- and outgoing 
to the singularity of the cardioid at 0j = ivr because 



9L(0i_i,0i) dL{(j)i,(f)^ 



+ 



i+l) 



. 



(5) 



bi=±TT 



In order to find cusp orbits we can therefore look for critical points of an action with fixed 
initial and final point 



WliXn] a, u) = L{a, xq) + L{xo, xi) H h cu) . 



(6) 
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Note that describes an orbit with n reflections not counting reflections at a and 
uj. It is sufficient to restrict to cusp orbits hitting the cusp just once, as these can be 
combined to multiple cusp orbits by concatenation. Similar to the case of regular orbits 
in the cardioid we can find cusp orbits by looking for maxima of W^. The proof is as for 
WP, the only difference is that (^o) is replaced by L(0„_i, u) + L{a, (po). For the 

cardioid with a = u = tt the Hessian of this function is diag(— cos0o, — cos (f)n-i) ■ Thus 
it is negative definite for orbits which hit the cusp from the right (both angles, before 
and after hitting the cusp, have modulus smaller than tt/2) and the same proof as in 



22| works. Unfortunately we have not been able to show that orbits with modulus of 0o 
or (pn-i larger then 7r/2 are also minima. Even if such an orbit is initially pruned, e.g. 
A^"'CC for n > 2, it can become physical under continuation. The shortest code for which 
we cannot find a corresponding minimum in the action is A^B^CC. However, using the 
symmetry line method described in to find this orbit shows that it is in fact s-pruned. 

The initial conditions for the numerical determination of the maximum of for a 
given cusp orbit are obtained similarly to the case of regular orbits: first the code word is 
rotated so that it starts with CC. The same procedure as above is applied for intermediate 
(f)i, except in the end we just put (f)^ = z = . To follow orbits under parameter variation 
we use a predictor-corrector method, see e.g. [^]. For a slightly changed parameter 
e + Ae an initial guess for the orbit is computed in the predictor step, which is used in 
the corrector step to find the actual orbit at the new parameter. We continue the orbit 
$ = (00, • • • , 0n-i) under parameter variation in e by expanding 

VW^P($ + A<l>; e + Ae) ^ WW^ + D^W^/^^ + . (7) 

Because $ is an orbit at parameter e we have VW^ = 0. Hence we can solve for A$ 
provided that D'^W^ is nonsingular. This gives the predictor $ + A$ for the continued 
orbit at parameter e + Ae. The corrector step consists in Newton's method taking the 
predictor as the initial guess for the orbit at the new parameter e + Ae:. Note that for 
continuing orbits away from the cardioid it is necessary to use Newton's method, because 
for e < 1 not all the critical points of are maxima. If the orbit cannot be found 
for the new parameter e + Ae the step size is decreased. If under parameter variation 
det D^W^ = is reached, we have found the point (in most cases, see Sec. ||) at which 
the orbit is created. In order to be able to continue orbits through this singular point in 
the parametrization by family parameter e we consider the augmented problem, see e.g. 
29| , chpt. 4.4], in which the parameter is itself a variable, and the continuation uses the 



arc length along the curve in ($, e) space as a new parameter. 

A cusp orbit of length n has n — 1 points ^n~i = (V'o, • • • , V'n-2) and one point at 
the cusp, ^n-i satisfies the equation VVF^_i(\l'n-i; tt, vr) = at £ = 1. To continue 
a cusp orbit away from the cardioid we make the fixed angle vr in PV^_]^(\l'„_i; tt, tt) an 
additional variable and extremise W^p($„) for e < 1, where initially = (^E'n-i, tt). Note 
that because of (^) we have Viy^($ri) = at £ = 1 as before. However, the Hessian 
of is not definite at this point. For the continuation the important fact is that it is 
nondegenerate, which we will now show. First note that 



d(f)id(j)j 



=±7r 



COS (pi^i + COS </)j+i if i = j 
else 



and that outside the cusp the second derivatives of and coincide. Therefore we 
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find 

('& \ 

) t9 = COS ?/'o + COS ?/'„_2 . 

We have already shown that D'^Wl_i is negative definite, and we also noted that for 
physical orbits the cosines in -(9 are both positive. Therefore D^W^ is nondegenerate at 
cusp orbits. The difference to orbits that do not hit the cusp is that here the Hessian is 
not negative definite, instead it has one positive eigenvalue. In any case these orbits can 
be continued away from the cardioid. By obvious extension of the above the case where 
a finite orbit hits the cusp more than once can be included as well. 

It is surprising that the continuation of cusp orbits away from the cardioid is not 
singular, because these orbits are "infinitely unstable" in the cardioid. This can be seen 
using a formula for the residue R of an orbit ||30[ 



IdetDWP , „ d^WP 
R= -- , where Bi - 



4 nr=o B^ 

with cyclic boundary conditions. Observe that for e — 1 we have that Bq 0, -Bn-i 
due to (H) and thus R diverges. The limiting behaviour can be found by expansion around 
e = 1. This gives Bi = {1 — e) cos (pi for i = and i = n — 1, so that we find 

1 detDW^_i / 1 



R ~ -TT^ ^——^r^ + 



4(1-^)' K=lB^ V COS 00 cos 

near 6 = 1. Note that for multiple cusp orbits a similar result holds with the power of 
1 — e decreased by one for each point in the cusp. 



3 Bifurcations of orbits 

To study the bifurcations of a periodic orbit we start in the cardioid and continue it from 
e = 1 towards the circle e = as long as it exists. We include all cusp orbits and also 
all pruned orbits, because they typically become physical for smaller e when the billiard 
becomes closer to being convex. An orbit can cease to exist if it collides with another 
one at residue R = 0, corresponding to an eigenvalue 1 of the linearised Poincare map. 
For all other finite R orbits cannot disappear, which can be seen by applying the implicit 
function theorem to (|^. This only holds for smooth maps, so the singular bifurcations 
which occur in a non-convex billiard when there is a tangency of an orbit segment and 
the boundary can occur with any residue. The variational approach is blind to these 
singular bifurcations, except if there is a refiection at the point of tangency. In this case 
the stability of the orbit diverges at the tangency. Numerically the continuation is done 
for decreasing e. However, when describing bifurcations we always consider an increasing 
parameter. The following possibilities arise: 

• The orbit exists down to e = 0. In the integrable circle billiard periodic orbits can 
be labeled by their rational rotation number u = p/q- Under small perturbations 
the corresponding rational torus breaks into a Poincare-Birkhoff chain consisting of 
one unstable and one stable orbit of period q. The existence of rotational p/q orbits 
for larger e is guaranteed by Aubry-Mather theory [^T], (see ||2^ for a review) 
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Figure 1: The trace TrM = —AR + 2 of the monodromy matrix is shown for the shortest 
nonsymmetric periodic orbit AAABABB and the cusp orbit AABABCC in the cardioid under 
variation of the parameter e. At e = 0.817422 . . . these periodic orbits are created in a saddle- 
node bifurcation (sn). Also shown are the corresponding periodic orbits at e = 0.85. 




Figure 2: Under parameter variation the AB orbit has a symmetry breaking period tripling 
bifurcation (sl/3) at e = 0.3728... creating the unstable orbit BACCCC and the initially 
elliptic orbit BAABCC. At e = \/2 — 1 a period doubling (pd) bifurcation gives birth to the 
AABB orbit. This orbit has full symmetry and therefore has a tangency with TrM = —2 at 
e = 0.4587 .... 
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Figure 3: Example of another typical bifurcation scenario for orbits with symmetry. First at 
e = 0.6953235 . . . the two orbits AAABB (stable) and ACCCC (unstable) are created in a 
saddle-node bifurcation. Notice that the orbit AAABB becomes inverse hyperbolic but again 
elliptic at e = 0.697988 .... At e = 0.709705 ... a symmetry breaking pitchfork bifurcation (pf) 
takes place in which the elliptic AABCC (and its symmetric partner BAACC) are created. All 
orbits are shown for e = 0.8. 




0.675 0.676 0.677 0.678 0.679 e 0.680 



Figure 4: The bifurcation diagram for the unstable pitchfork bifurcation (upf) is the mirror 
image of the previous one in Fig. ^. Here either partner in the saddle-node bifurcation {e = 
0.675862 . . . ) eventually becomes stable. The initially unstable orbit BAAABCC becomes stable 
at £ = 0.678499 . . . and creates BAACCCC and its symmetric partner AABCCCC (not shown). 
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as long as the billiard is strictly convex, i.e. for e < 1/2. The unstable p/q orbits 
maximise the action W, the (initially, i.e. for e > 0, e small) stable ones are called 
minimax, they are saddles of W with one up direction. The minimax orbits up to 
period 7 together with the e value at which they turn unstable are listed in Table |l|. 

The orbit only exists down to a nonzero parameter value at which it has R = and 
dR/de 7^ 0. In systems with a reversing symmetry generically one of the following 
bifurcations occurs; the orbit collides with 

1. another orbit with opposite sign of i? in a saddle node bifurcation (sn); there 
is no parent orbit (Fig. 1). 

2. another orbit with the same sign of i? in a symmetry breaking pitchfork bi- 
furcation (pf); the parent orbit has i? = (Fig. 3). The subcritical pitchfork 
bifurcation in which the parent turns stable is denoted by upf (Fig. 4). 

3. itself in a period doubling bifurcation (pd); the parent orbit has R = 1. The 
subcritical period doubling in which the parent turns stable is denoted by upd 
(Fig. 2). 

4. another orbit with opposite sign of i? in a p/g-tupling bifurcation, g > 3; 
the parent orbit has R = sin^^np/q). The "touch and go" bifurcation with 
dR/de = for q = 3 and possibly 4 does not terminate the orbit because it 
also exists for smaller e. 



The above cases exhaust the list of generic bifurcations with one symmetry |3^. The 
family of limagon billiards has two reversing symmetries, X and T, see [Q. Therefore 
there are more generic bifurcations of symmetric orbits in this family p^ , p5| , |36| . E.g., 
the period tripling of an orbit with full symmetry creates four orbits with less symmetry 
whose phase portrait looks like a six-tupling bifurcation. We denote this bifurcation by 
sl/3 (Fig. 2). 

All bifurcations for periodic orbits up to period 7 are listed in Table 0. The parent orbit 
(which does not exist for sn) exists before and after the bifurcation. With the exception 
of the g-tupling bifurcations it changes its stability at the bifurcation, i.e. it turns from 
elliptic to hyperbolic (for pf, pd) or vice versa (upf, upd). If there are two partners in a 
bifurcation the stable one is listed first. In the last column of Table |^ the width of the 



p/q 


code 


e 


S-tangency 


1/2 


CC 


0.41421356237 




1/3 


CAC 


0.36810610044 




1/4 


CAAC 


0.35203466109 




1/5 


CAAAC 


0.34773961433 




2/5 


CABAC 


0.30229901387 




1/6 


CAAAAC 


0.34830176077 


0.89431500512 


1/7 


CAAAAAC 


0.35095913185 


0.81112873450 


2/7 


CAABAAC 


0.29265846151 




3/7 


CABABAC 


0.27009692367 





Table 1: Minimax orbits with winding number p/q, the corresponding code in the cardioid 
billiard and the e value for which they have R = 1 and turn unstable. The last column lists 
the e value where the orbit has a tangency with the boundary and becomes unphysical 
for larger e. 
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parent 


type 


sym 


partner 


partner 


e 


Ae 


CC 


pd 


full 


AABB 




0.41421356237 


4.45e-2 






full 


AABB 




0.45873253228 


7.34e-2 


AABB 


pf 


T 


CABC 


CBAC 


0.53208888624 


2.30e-2 




sn 


X 


AAABB 


ACCCC 


0.69532356899 


1.87e-3 






X 


AAABB 




0.69798858052 


1.17e-2 


AAABB 


pf 




CAABC 


CBAAC 


0.70970577826 


4.58e-3 




sn 


X 


AAAABB 


AACCCC 


0.79431559532 


5.55e-4 






X 


AAAABB 




0.79572052126 


3.33e-3 


AAAABB 


pf 




AAACCB 


AAABCC 


0.79904759065 


1.43e-3 


CAC=l/3 


pd 


X 


AAABAB 




0.36810610044 


2.84e-2 






X 


AAABAB 




0.40610662358 


3.92e-2 


AAABAB 


pf 




ABAACC 


AABACC 


0.44533792566 


1.21e-2 




sn 


full 


AAABBB 


ACCBCC 


0.86753591308 


3.61e-4 






full 


AAABBB 




0.86789693974 


1.41e-3 


AAABBB 


pf 


T 


BBAACC 


AABBCC 


0.86931082693 


5.39e-4 




sn 


T 


AABABB 


ABABCC 


0.66948391241 


1.79e-4 


CC 


sl/3 


T 


BACCCC 


ABCCCC 


0.37282469327 





CC 


sl/3 


X 


BAABCC 


ABBACC 


0.37282469327 


6.12e-2 




sn 


X 


AAAAABB 


AAACCCC 


0.84879163221 


2.31e-4 






X 


AAAAABB 




0.84943955625 


1.35e-3 


AAAAABB 


pf 




BAAAACC 


AAAABCC 


0.85079348228 


5.74e-4 




sn 


X 


AAAABAB 


CAACCAC 


0.55897639858 


2.10e-3 






X 


AAAABAB 




0.56209428854 


8.98e-3 


AAAABAB 


pf 




ABAAACC 


AAABACC 


0.57107790543 


3.44e-3 




sn 


X 


AAAABBB 


CAACCBC 


0.93351335014 


6.67e-5 






X 


AAAABBB 




0.93359478150 


2.57e-4 


AAAABBB 


pf 




BBAAACC 


AAABBCC 


0.93385152715 


9.61e-5 




sn 


X 


AABAABB 


CABACCC 


0.61088771114 


1.83e-4 






X 


AABAABB 




0.61991127863 


3.88e-3 


AABAABB 


pf 




BAABACC 


ABAABCC 


0.62379601274 


1.86e-3 




sn 


X 


ACCCCCC 


BAAABCC 


0.67586228774 


7.77e-6 






X 


BAAABCC 




0.67887512963 


-3.76e-4 


BAAABCC 


upf 




BAACCCC 


AABCCCC 


0.67849936105 







sn 




AAABABB 


AABABCC 


0.81742214822 


1.46e-5 




sn 




BACCACC 


AABBACC 


0.75730213644 


2.35e-7 



Table 2: Table of bifurcations for orbits up to period 7 in the family of limacon billiards. 
Since maximizing orbits do not bifurcate they are not included. Note the sequences of 
symmetric saddle nodes followed by a pitch fork bifurcation, see Fig. 3,4. At the end of 
the table there are the first two pairs of orbits born in a saddle node bifurcation without 
symmetry, see Fig. 1. For higher periods most orbits are created in this type of bifurcation. 
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stable interval is listed, where a negative number means that it extends to the left from 
the listed value of e. For symmetric orbits usually one of the partners has a second interval 
of stability. These orbits are listed as a single partner without a bifurcation type together 
with the value e for which R = 1. The same families will also be parents of pitch-fork 
bifurcations when i? = 0. This is listed next, so that for most symmetric orbits there 
are three consecutive entries that describe the essentials of the corresponding bifurcation 
tree (Fig. 3,4). For orbits with full symmetry the two intervals of stability join at a point 
where the graph of R{s) has a tangency with 1 (Fig. 2). In this way all stable intervals 
of orbits up to period 7 and their bifurcations are listed in Table ^. 

The relation between the winding number p/q and the symbolic code for the maxi- 
mizing and minimax orbits is simple. The code of the maximizing p/q orbit is given by 
the word with symbols 




A if [ip/q] = [{i + l)p/q] 
B otherwise 



where [x] denotes the integer part of x and i = 0, . . . ,q — 1. The corresponding minimax 
orbit of the same rotation number is obtained by replacing the first and the last letter by 
C, see Table |I] for some examples. With these rules we can predict at e = 1 which orbit 
is going to have a particular rotation number at the integrable limit. 

A typical simple saddle node bifurcation is shown in Fig. |I[ Depending on the symme- 
try we observe more complicated sequences of bifurcations, see Fig. 0, Fig. ||, and Fig. | 
for examples. These are the four basic types of (sequences of) bifurcations that occur up 
to period 7, compare to the corresponding entries in Table |^. 



4 Statistics of Bifurcations 

We numerically computed all stable intervals of all orbits up to period 13. This is a 
complicated task that involves continuing 3697 periodic orbits. In particular it came as a 
surprise that the width of the stability intervals for some of the orbits could not be calcu- 
lated using double precision, i.e. approximately 15 significant digits. E.g. already at pe- 
riod 10 the smallest interval of stability occurs for the saddle node pair AAACCBCCCC, 
AAAABBCCCC , and is only Ae ^ 4.1 ■ 10"^*^. In order to be able to calculate such small 
intervals we used the doubledouble package p^, which offers 30 significant digits. This 



enabled us to determine all stable intervals of continued orbits up to period 13. Studying 
higher period would require even more precision. We now report on the results of our 
numerical study. 

Up to period 7 the results are listed in Table 2. Fig. ^ shows the width Ae of intervals 
of stability as a function of e (at the left boundary of the interval) for orbits of period 11 
and for all orbits up to period 13. The number of points grows with the period like the 
number of orbits (including cusp orbits) in the cardioid, i.e. exponentially with topological 
entropy ln(l + \/2) |2^. Smaller intervals occur for larger e. The trivial upper bound 
is ^s{e) = 1 — e. For fixed period a lower bound is approximately given by a straight 
line. There are conspicuous patterns that can be related to the structure of the code 
word. The detailed structure is quite complicated, and it is difficult to decide whether 
for each parameter there exist elliptic orbits. Even though the number of stable intervals 
grows exponentially with the period this does not mean that they are easily visible in 
phase space. In [O, O it is referred to an unpublished numerical study by Li and Robnik 



11 







10 
Ae 



10 



10 



10 



10 



10 



15 



■20 



■25 



- 


+• * 


, + 

+ +++ * 

+ + • + ■* + + + . 

4. + + ' • i +• + • 


- 




. V*. \* a *" ; 

• • • J. - 


- 




• *. , 

•••^^-. v : • 

• • • ' ': ' 
•* • • 


- 

: 


n=ll 


• • - 
• • 




symmetric 
not symmetric 
l_e 


• 




1 1 


1 1 







0.2 



0.4 



0.6 



0.8 




Figure 5: Width Ae of stable intervals as a function of e for orbits of period 11 and for 
orbits up to period 13. 
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Figure 6: Averaged widths of stability intervals for symmetric and non-symmetric orbits up to 
period 13. The slopes are approximately equal to the negative of the topological entropy. 
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Figure 9: The plot of the parameter of creation of A^BB orbits as a function of n shows 
that there are elliptic orbits created in saddle node bifurcations arbitrarily close to the cardioid 
billiard. 
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(1994), giving that if there is a stabihty island for e > 0.555 its relative area has to be 
smaller than 5 • 10^''. This is in agreement with examples considered in 0. 

To obtain a quantitative estimate of the behaviour of width of stability intervals with 
the period we now consider the average Ae of Ae for each period, see Fig. |^. Here 
it is important to treat orbits with and without symmetry separately because orbits 
with symmetry tend to have larger intervals of stability. A reason for this is that many 
symmetric orbits undergo symmetry breaking pitchfork bifurcations, see Fig. ^. Therefore 
they have to intersect TrM = 2 with finite slope, so that the second interval of stability is 
larger than the first. There is a difference between even and odd periods. On average non- 
symmetric orbits with odd period have smaller stable width than those with even period. 
This is mainly due to the fact that there are no period doubled orbits among them, which 
tend to have a larger stable width because of their finite slope in the bifurcation diagram 
at creation (see e.g. Fig. The average width of the parameter interval on which the 
orbits are stable follows a power law with the period. For the non-symmetric orbits the 
exponent is close to the topological entropy 0.881. The actual growth of the number of 
intervals of stability is shown in Fig. ^ Therefore the width of the intervals decreases at 
the same rate as the number of orbits increases, such that it is not possible from these 
estimates to determine the asymptotic behaviour of the sum of all intervals. For the 
symmetric orbits a approximately half the exponent is expected and measured. 

We can conclude that the behaviour of this family of billiards is much more complicated 
than one would expect on first sight. Even arbitrarily close to the ergodic limit of the 
cardioid there exist tiny elliptic islands. This can be seen by considering the family of 
orbits A^BB that limits onto the whispering gallery orbit. This family illustrates the non- 
uniform hyperbolicity of the billiard map, its Liapunov exponent tends to as tends to 
infinity, see . In Fig. ^ a plot of the parameter En where A^BB is born in a saddle node 
bifurcation is shown. The data is well described by the fitted curve e = 1 — arr'^l'^ where 
a = 1.65. This shows that arbitrarily close to the cardioid there are stable elliptic orbits 
created in saddle node bifurcations. In Fig. ^ this family of orbits appears approximately 
as an upper bound for Ae, i.e. they tend to have the largest stability intervals for given e. 

Whether there are parameter values for which the corresponding limagon billiard is 
ergodic cannot be decided by a numerical experiment like the one we performed. The 
results indicate, however, that the set of parameter values for which the map is ergodic 
(if larger than one point) has a very complicated structure. 

To our knowledge this is the first time that the scaling of stability intervals with the 
period has been extensively studied. It is well known that in hyperbolic systems there 
is a scaling of the instability of periodic orbits with the period. This expresses the fact 
the the phase space volume is a constant so that growth of instability has to be exactly 
compensated by growth of the number of orbits |3^. Our results indicate that there is a 
similar relation for stability intervals in parameter space for families limiting to hyperbolic 
system. 
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